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driver  wave  modeling.  The  analysis  is  based  on  the  concept  that  river  wave  behavior  is  determined  by  the  balance  be¬ 
tween  friction  and  inertia.  The  Saint-Venant  equations  are  combined  to  form  a  system  equation  that  is  written  in 
dimensionless  form.  The  dominant  terms  of  the  system  equation  change  with  the  relative  magnitudes  of  a  group  of 
dimensionless  scaling  parameters  that  quantify  the  friction/inertia  balance.  These  scaling  parameters  are  continuous, 
indicating  that  the  various  types  of  river  waves  and  the  transitions  between  them  form  a  spectrum.  Additional  data 
describing  the  physical  variability  of  rivers  and  waves  are  incorporated  into  the  analysis  by  interpreting  the  scaling 
parameters  as  random  variables.  yThis  probabilistic  interpretation  provides  an  improved  estimate  of  the  friction/inertia 
balance,  further  insight  into  the  continuous  nature  of  wave  transitions,  and  a  measure  of  the  reliability  of  wave  type 
assessments  near  a  transition.  Case  studies  are  used  to  define  the  scaling  parameter  ranges  that  represent  each  wave 
type  and  transition,  and  to  provide  data  with  which  to  evaluate  the  usefulness  of  the  analysis  for  general  application,  i 
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NOMENCLATURE 

P\y<x] 

probability  that  random  variable  y  is 

less  than  x 
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hydraulic  radius  of  the  channel 

^0 

initial  wave  amplitude 
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B 
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St 
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c 
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flow 

cd 

celerity  of  diffusion  waves 

S0 

river  bed  slope 

celerity  of  kinematic  waves 

Sx 

standard  deviation  of  random  variable  x 

'-m 

measured  wave  celerity 

t 

time 

C, 

Courant  number,  scaling  parameter  of 

t* 

dimensionless  time 

the  Saint-Venant  system  equation 

V 

flow  velocity 

c. 

dimensionless  Chezy  conveyance  co¬ 

vo 

velocity  scaling  variable 

efficient 

V* 

dimensionless  velocity 

D 

dimensionless  diffusion  coefficient 

K 

coefficient  of  variation  of  random  vari¬ 

D, 

scaling  parameter  of  the  Saint-Venant 

X 

able  X 

system  equation 

X 

longitudinal  distance 

E\y ] 

expected  value  of  random  variable 

x • 

dimensionless  distance 

fly) 

0  distribution  of  random  variable  y 

y 

flow  depth 

Fi,Fl 

friction  scaling  parameters  of  the 

y0 

depth  scaling  variable 

Saint-Venant  system  equation 

y* 

dimensionless  depth 

F0 

Froude  number 

parameters  of  a  0  distribution 

g 

acceleration  due  to  gravity 

A t 

time  scaling  variable 

i 

V-T 

Ax 

length  scaling  variable 

k 

variable  relating  depth  to  hydraulic 

1 

coefficient  of  exponential  amplitude 

radius;  arbitrary  constant 

decay  of  reservoir  waves 

kltk2 

constants 

P 

correlation  coefficient 

ANALYSIS  OF  RIVER  WAVE  TYPES 


Michael  G.  Ferrick 


INTRODUCTION 

Long-period  waves  in  rivers  are  a  consequence  of 
unsteady  flow,  which  may  occur  as  a  result  of  hydro¬ 
electric  power  generation  or  flow  control  at  a  dam, 
the  breach  of  a  dam,  the  formation  or  release  of  an 
ice  jam,  or  rainfall/runoff  processes.  River  waves 
have  flow  depths  that  are  several  orders  of  magni¬ 
tude  smaller  than  their  wavelengths,  so  they  are 
classified  as  shallow-water  waves.  The  study  of  river 
waves  is  complicated  by  the  hydraulic  properties  of 
rivers  and  the  time  scales  of  the  wave  motion,  which 
range  over  several  orders  of  magnitude.  Natural  riv¬ 
ers  vary  greatly  in  size  and  in  discharge.  Impounded 
rivers  are  both  wider  and  deeper  than  they  were  in 
their  natural  state,  and  they  have  diverse  hydraulic 
properties.  Flood  waves  in  large  river  systems  de¬ 
velop  over  several  days,  while  those  in  some  small 
rivers  can  form  in  several  minutes.  Flow  control  and 
hydroelectric  turbine  gates  may  be  adjusted  rapidly, 
creating  abrupt  flow  waves  both  upstream  in  the 
reservoir  and  downstream  in  the  tailwater  river.  Be¬ 
cause  of  this  diversity,  general  physical  insights  re¬ 
garding  river  waves  are  elusive. 

Observations  indicate  that  river  waves  are  a  pri¬ 
mary  cause  of  river  ice  cover  breakup.  Therefore, 
a  theory  describing  the  interactions  between  a  river 
wave  and  the  ice  cover  would  have  great  practical 
importance.  When  a  river  becomes  covered  with  ice, 
frictional  energy  dissipation  increases,  but  the  effect 
of  this  increased  friction  on  river  wave  behavior  has 
not  been  established.  This  understanding  is  a  pre¬ 
requisite  to  the  development  of  theories  of  ice  cover 
stability  in  the  presence  of  river  waves. 

The  Saint-Venant  equations  of  continuity  and 
momentum  are  generally  used  to  describe  unsteady 


flow  in  rivers.  Within  this  basic  model,  river  waves 
may  be  classified  as  dynamic,  gravity,  diffusion,  or 
kinematic  waves,  corresponding  to  different  forms 
of  the  momentum  equation.  Dynamic  waves  are  de¬ 
scribed  by  retaining  all  terms  of  the  momentum  equa¬ 
tion.  As  is  typical  of  shallow-water  waves,  dynamic 
river  waves  have  a  celerity  that  is  related  to  the  water 
depth.  The  momentum  equation  for  gravity  waves 
ignores  the  effects  of  bed  slope  and  viscous  energy 
losses;  gravity  waves  propagate  at  the  dynamic  wave 
celerity  and  their  flows  are  dominated  by  inertia. 
Diffusion  and  kinematic  waves  are  at  the  opposite 
extreme.  The  diffusion  wave  momentum  equation 
ignores  inertia,  and  the  kinematic  wave  equation  ig¬ 
nores  both  inertia  and  the  pressure  gradient  caused 
by  varying  flow  depths  over  distance.  Both  of  these 
wave  types  travel  at  a  celerity  governed  by  the  fric¬ 
tional  resistance  of  the  river  bed.  This  celerity  is  re¬ 
lated  to  the  velocity  of  the  flow  and  is  significantly 
slower  than  the  dynamic  wave  celerity.  Diffusion 
and  kinematic  wave  propagation  requires  the  move¬ 
ment  of  a  large  quantity  of  water.  In  this  paper 
these  waves  are  frequently  grouped  together  and  are 
called  bulk  waves. 

Our  present  understanding  of  the  relationships 
between  various  river  wave  types  is  based  on  linear 
stability  theory  (Ponce  and  Simons  1977,  Menendez 
and  Norscini  1 982).  This  linear  theory  provides  use¬ 
ful,  but  largely  qualitative,  insights  into  the  behavior 
of  each  wave  type,  and  the  transitions  between  wave 
types  are  not  considered:  waves  of  a  given  type  in¬ 
clude  or  exclude  specific  terms  or  processes  in  the 
momentum  balance.  The  concept  of  step  changes 
between  river  wave  types  is  not  reasonable.  There 


must  be  cases  that  are  intermediate  between  wave 
types,  suggesting  the  need  for  a  treatment  of  wave 
transitions.  Wave  types  are  identified  either  by  us¬ 
ing  the  linear  theory  or  by  comparing  magnitudes 
of  normalized  terms  of  the  dynamic  wave  momen¬ 
tum  equation  (Henderson  1963,  Woolhiser  and 
Liggett  1967),  and  then  applying  judgment  to  in¬ 
terpret  the  results. 

The  object  of  this  paper  is  a)  to  develop  a  quanti¬ 
tative  method  for  identifying  river  wave  types  and 
b)  to  clarify  the  relationships  between  wave  types. 
The  analysis  is  based  on  the  principle  that  the  bal¬ 
ance  between  friction  and  inertia  determines  river 
wave  behavior.  The  Saint-Venant  equations  are 
combined  to  form  a  system  equation.  Written  in 
dimensionless  form,  the  system  equation  provides 
scaling  parameters  that  quantify  the  magnitudes  of 
all  terms  in  the  equation  and  indicate  the  relative 
importance  of  friction,  inertia,  and  pressure  gradient 
effects  on  a  wave.  By  interpreting  the  scaling  param¬ 
eters  probabilistically,  additional  data  can  be  incorpo¬ 
rated  into  the  analysis  to  account  for  the  variable 
physical  characteristics  of  a  river  and  a  wave.  This 
more  general  interpretation  provides  an  improved 
estimate  of  the  friction/inertia  balance,  insight  into 
the  continuous  nature  of  transitions  between  wave 
types,  and  a  measure  of  the  reliability  of  wave  type 
assessments  near  a  transition.  Finally,  we  identify 
the  scaling  parameter  ranges  that  correspond  to  each 
wave  type  and  transition  with  data  from  case  studies. 
These  case  studies  encompass  a  wide  range  of  river 
and  wave  conditions,  attesting  to  the  general  utility 
of  the  approach. 


BACKGROUND 

C'ourant  and  Friedrichs  (19  '  ( noted  the  analogy 
between  nonlinear  wave  motion  in  gases  and  in  shal¬ 
low  water,  and  they  developed  the  mathematics  to 
treat  fust-oidei .  quaai-linear  hyperbolic  flow  equa¬ 
tions  for  funt'.ona  ot  two  independent  variables. 

( ontmuity  ot  all  (unctions  and  all  required  deriva¬ 
tives  ol  these  lumtiom  was  assumed  Stoker  ( 1957) 
applied  thn  theorv  f«>  the  studv  of  river  flow  waves, 
pro*whng  the  mathematical  and  omceptuaJ  baas  for 
oui  preseni  under sunding  *  nvet  waves.  Following 
Stidtet  mt  .onssdr  idealised  rivet  with  a  wide 
retUmptUi  pstsenatn  ixanori  ud  no  local  inflow 
The  Sunt  Vena n i  rvjuai cmi  ot  mass  conservation 
and  momentum  tvalar.r  tie  then  written 


dr 


dx 


(2) 


where  y 
v 


g 

So 

Sf 

x 

t 


the  flow  depth  (m) 
the  flow  velocity  (m/s) 
acceleration  due  to  gravity  (m/s2 ) 
the  river  bed  slope 
the  energy  gradient  of  the  flow 
the  longitudinal  distance  (m) 
time  (s). 


Stoker  transformed  eq  1  and  2  into  their  character¬ 
istic  form: 


v  +  2e+g(S{-S0)t=kl  (3) 

along  =  v  +  c. 


and 


v  -  2c  +g(Sf  -  S0)  t  =k2  (4) 


i  dx 

along  —  =  v -c, 
at 


where  c  -  yjgy  defines  the  shallow-water  surface 
wave  celerity  and  kt  and  k2  are  constants.  The  solu¬ 
tions  of  dx/dt  -v±c  yield  two  distinct  sets  of  curves 
in  the  x-t  plane,  called  characteristics.  The  inverse 
slopes  of  these  curves  define  the  dynamic  wave  celer¬ 
ity,  that  is,  the  speed  of  shallow-water  waves  in  a 
flowing  stream.  The  dynamic  wave  celerity  is  inde¬ 
pendent  of  the  river  bed  and  energy  slopes,  and  ap¬ 
pears  to  be  the  only  wave  speed  contained  in  the 
governing  equations. 

The  momentum  equation  describing  a  gravity 
wave  is  obtained  from  eq  2  by  omitting  the  river 
bed  slope  and  energy  gradient  terms.  The  theory  of 
characteristics  can  also  be  applied  for  these  waves, 
and  the  resulting  equations  are  simplified  forms  of 


eq  3  and  4: 

v  +  2c  =  k, 

along 

dx 

Si-”*'- 

(5) 

v-2 c  =  *, 

along 

dx 

dt  =v~c’ 

(6) 

6*  ‘U 


(1) 


<i\ 
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Waves  described  by  eq  5  and  6  have  been  termed 
“ample”  waves.  Gravity  waves  are  undamped  and 
propagate  at  the  dynamic  wave  celerity. 


The  diffusion  wave  momentum  equation  is  ob¬ 
tained  by  neglecting  the  inertia  terms  [3v/3r,  v(3v/ 
dor)]  in  eq  2.  Differentiating  both  this  momentum 
equation  and  eq  1  with  respect  to  x,  and  differen¬ 
tiating  the  momentum  equation  with  respect  to  t, 
yields  three  equations  that  are  combined  with  eq  1 
to  eliminate  derivatives  of  the  depth.  The  Chezy 
equation,  with  dimensionless  conveyance  coefficient 
C,  assumed  constant,  is  used  to  describe  the  energy 
slope.  Performing  these  operations  yields 


3v^  dv  nd2v 

a7+c<*ai  =  Da? 


(7) 


where 


gCi  RS0 


c^V- 


and  D- 


gCi  Ry 
2v  • 


Equation  7  is  an  advective-diffusion  equation  with 
R  the  hydraulic  radius  of  the  channel.  If  we  retain 
the  wide-channel  assumption,  the  hydraulic  radius 
is  equivalent  toy  Ik,  where  k  =  1  for  open  water  con¬ 
ditions  and  k  =  2  if  the  channel  is  covered  with  ice. 

The  kinematic  wave  momentum  equation  is  ob¬ 
tained  by  omitting  the  inertia  and  the  derivative  of 
depth  with  distance  (dy/dx)  terms  in  eq  2.  Then, 
proceeding  as  with  diffusion  waves  yields  the  well- 
known  kinematic  wave  equation 


3v  3v  . 
3/  +  C|t  3x  =  0  ’ 


(8) 


where 


ck  =  2  v  ’ 


The  important  distinction  between  the  bulk  waves 
is  that  diffusion  waves  attenuate,  while  kinematic 
waves  do  not. 

Comparing  eq  3  through  8,  we  find  that  the  equa¬ 
tion  for  river  wave  celerity  changes  with  the  form  of 
the  momentum  equation  that  is  used.  The  celerities 
of  diffusion  and  kinematic  waves  are  both  related  to 
the  velocity  of  the  flow.  The  celerities  of  these  bulk 
waves  are  typically  much  smaller  than  the  dynamic 
wave  celerity  and  are  generally  lumped  together  and 
termed  the  kinematic  wave  celerity. 

Stoker  (19S7)  argued  that  because  the  equations 
describing  dynamic  waves  are  more  general,  imper¬ 
ceptible  dynamic  waves,  or  “forerunners,”  must 
occur  simultaneously  when  the  primary  wave  is  kine¬ 


matic.  This  concept  of  an  ever-present  role  for  dy¬ 
namic  waves  in  river  flow  mechanics  is  widely  accept¬ 
ed.  Therefore,  river  wave  studies  based  upon  “sim¬ 
plified”  gravity,  diffusion,  and  kinematic  wave  equa¬ 
tions  are  perceived  as  inherently  less  accurate  than 
those  using  the  more  complete  dynamic  wave  equa¬ 
tions. 

The  governing  equations  for  each  wave  type  are 
nonlinear,  and  the  geometry  of  natural  channels  is 
highly  variable.  These  conditions  generally  require 
that  numerical  solutions  be  used  to  describe  rive, 
waves.  The  algorithms  typically  chosen  to  solve  the 
dynamic  wave  equations  in  applications  use  first-  or 
second-order  accurate  difference  approximations 
and  coarse  numerical  meshes.  The  solutions  obtained 
with  these  models  therefore  include  significant  errors, 
due  to  truncation  and  discretization,  that  cause  nu¬ 
merical  diffusion  and  dispersion  of  the  solution. 

When  the  effects  of  certain  physical  processes  are 
small  relative  to  other  processes,  the  terms  of  the 
governing  equation  have  widely  different  magnitudes. 
Equations  of  this  type  are  termed  “stiff,”  and  are 
difficult  to  solve  numerically.  Woolhiser  and  Liggett 
(1967)  reported  numerical  difficulties  that  resulted 
from  modeling  predominantly  kinematic  waves  using 
dynamic  wave  characteristic  equations.  The  stability 
of  explicit  numerical  methods  depends  on  the  dy¬ 
namic  wave  celerity,  even  for  cases  where  the  inertia 
terms  are  negligible.  When  stability  problems  occur, 
the  cure  is  often  achieved  by  increasing  numerical 
diffusion  in  the  algorithm,  further  degrading  the  ac¬ 
curacy  of  the  solution.  Taken  together,  these  con¬ 
siderations  imply  that  more  complete  equations  may 
not  yield  more  accurate  river  wave  simulations. 

The  problems  associated  with  universal  application 
of  the  dynamic  wave  equations  suggest  the  alternative 
approach  of  using  physical  insight  to  identify  an  ap¬ 
propriate  wave  type.  However,  wave  types  cannot  at 
present  be  identified  with  a  known  degree  of  certain¬ 
ty.  For  example,  Cunge  et  al.  (1980)  stated  that 
rapidly  varying  river  flows  require  the  use  of  the  dy¬ 
namic  wave  equations,  but  Ferrick  et  al.  (1984) 
studied  a  laije  number  of  instantaneous  flow  releases 
in  two  rivers  and  found  that  in  all  cases  the  inertia  of 
the  flow  was  negligible.  The  logic  linking  rapid  flow 
variation  with  inertia  is  clear  but  insufficient  to  en¬ 
sure  dynamic  wave  behavior.  The  converse  example 
is  that  slowly  developing,  long-period  floods  behave 
as  kinematic  waves.  A  more  quantitative  representa¬ 
tion  of  the  roles  of  friction  and  inertia  would  make 
it  possible  to  characterize  and  identify  river  wave 
types. 


ANALYSIS 


terms  are  due  to  inertia.  This  system  equation  has 
the  form 


River  waves  propagating  at  the  dynamic  and  kine¬ 
matic  wave  celerities  are  generally  observed  in  differ¬ 
ent  situations.  In  impounded  rivers,  where  depth  is 
much  greater  than  wave  amplitude,  observed  wave 
celerities  are  those  of  dynamic  waves  (Ferrick  and 
Waldrop  1977,  Ferrick  1979),  but  wave  propagation 
observed  in  free-flowing  rivers  is  more  commonly  at 
the  kinematic  wave  celerity  (Ferrick  et  al.  1984). 
These  observations  suggest  that  it  may  be  possible 
to  distinguish  between  bulk  and  dynamic  waves 
simply  by  measuring  and  classifying  wave  celerity. 

We  might  also  reason  that  river  waves  are  shallow- 
water  waves  of  extremely  long  wavelength.  A  depth 
parameter  nondimensionalized  by  wavelength,  which 
follows  from  the  simple  harmonic  small-amplitude 
wave  theory  of  ocean  gravity  waves,  differentiates 
between  deep-water  and  shallow-water  conditions. 
Extending  these  ocean  wave  classification  ideas  to 
rivers,  in  light  of  typical  wavelength  differences  be¬ 
tween  rapidly  and  gradually  varying  flows,  suggests 
that  a  dimensionless  depth  or  a  dimensionless  wave¬ 
length  may  be  indicative  of  wave  type. 

A  connection  between  the  form  of  the  momen¬ 
tum  equation  and  the  physical  characteristics  of  a 
river  and  wave  would  provide  a  means  to  identify 
wave  type  and  would  clarify  the  relationships  be¬ 
tween  waves.  Each  of  the  river-wave  types  that  have 
been  identified  is  described  by  equations  that  are 
subsets  of  the  Saint-Venant  equations.  To  consider 
the  behavior  of  solutions  of  this  equation  system,  we 
will  assume  that  the  functions  describing  depth  and 
velocity  possess  continuous  second  derivatives  and 
that  the  river  bed  slope  is  constant.  Then,  differen¬ 
tiating  eq  1  and  2  with  respect  to  distance  x,  differ¬ 
entiating  eq  2  with  respect  to  time  t,  and  combining 
these  three  equations  with  eq  1  and  2  to  eliminate 
the  depth  derivative  terms  yields  a  system  equation: 
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The  second-order  terms  in  second-order  differential 
equations  are  generally  of  principal  significance.  As 
B 2  -  AC  is  greater  than  zero,  eq  9  is  hyperbolic  with 
a  pair  of  characteristic  curves  (Hildebrand  1962)  de¬ 
fined  by 

A$  -  2B(jt)  +  C  =  °.  (11) 

Solving  eq  1 1  for  dxldt  again  yields  the  dynamic 
wave  celerity  v  ±  c  as  the  inverse  slope  of  the  charac¬ 
teristics. 

The  assumption  implicit  in  the  development  of 
eq  3  and  4  is  that  all  the  processes  represented  in  eq 
2  are  of  comparable  magnitude.  However,  from  the 
perspective  of  the  system  equation,  the  behavic  -i 
a  river  wave  must  depend  on  the  relative  contribu¬ 
tions  of  inertia,  friction,  and  pressure.  When  inertia 
is  significant,  dynamic  waves  are  of  primary  impor¬ 
tance  in  river  flow  mechanics,  as  indicated  by  eq  3 
and  4.  The  role  of  dynamic  waves  diminishes  with 
the  relative  importance  of  inertia.  Bulk  wave  be¬ 
havior  indicates  that  the  second-order  terms  of  eq  9 
are  dominated  by  the  quasi-linear  first-order  terms; 
that  is,  friction  dominates  inertia. 

An  evaluation  of  the  relative  magnitudes  of  the 
terms  in  eq  9  and  a  basis  for  interpreting  these  re¬ 
sults  would  provide  a  quantitative  measure  of  the 
importance  of  each  process  in  the  momentum  bal¬ 
ance  and  would  indicate  the  wave  type.  Writing  eq 
9  in  dimensionless  form  permits  an  assessment  of 
these  relative  magnitudes.  Introducing  parameters 
v0,y0,  Ax,  and  Af,  which  are  in  some  sense  charac¬ 
teristic  of  the  flow  and  wave  motion,  we  rewrite  eq 
9  in  terms  of  dimensionless  variables  v*  =  v/v0,y*  = 
y/y0,x*  =x/Ax,  and  t*  =  t/At  as 
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in  which  the  Chezy  equation  with  constant  C*  is 
again  used  to  describe  the  energy  slope.  The  second- 
order  terms  in  eq  9  follow  from  the  inertia  and  pres¬ 
sure  gradient  terms  in  eq  2.  The  quasi-linear  first- 
order  terms  result  from  the  energy  or  friction  slope 
and  the  bed  slope  terms,  and  the  nonlinear  first-order 
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To  proceed,  we  must  define  the  physical  scaling 
variables  v0,y0  Ax,  and  A/  and  evaluate  SQ  and  C, . 
Mean  velocity  and  depth  characterize  the  flow  and 
provide  velocity  and  depth  scales  v0  andy0.  The 
length  scale  that  is  important  in  the  development  of 
river  waves  is  related  to  the  wavelength.  With  mean 
flow  depth  as  the  depth  scale,  the  comparable  length 
scaling  parameter  Ax  is  the  half-wavelength,  and 
kA x/y0  appearing  in  parameters  Fl  and  Fc  is  a  di¬ 
mensionless  wavelength.  The  effect  of  friction  on  a 
wave  is  cumulative  over  the  propagation  distance. 

At  distances  less  than  a  wavelength  from  the  point 
of  wave  origin,  Ax  is  taken  as  half  the  wave  propaga¬ 
tion  distance.  At  is  the  time  required  for  the  wave 
to  travel  distance  Ax.  Therefore,  Ax/ A t  is  the  mea¬ 
sured  wave  celerity  cm .  These  three  variables  are 
evaluated  with  information  from  a  pair  of  gauging 
stations.  To  obtain  measured  wave  celerity,  the  river 
distance  between  the  gauges  is  divided  by  the  elapsed 
time  between  the  wave  arrivals.  Wavelength  is  then 
obtained  by  multiplying  wave  celerity  by  the  mean 
time  for  wave  passage  at  the  gauges.  As  the  Chezy 
conveyance  coefficient  is  depth-dependent,  the  most 
representative  value  of  C.  for  a  reach  is  obtained 
from  steady-flow  stage  measurements  near  the  mean 
of  the  expected  range.  Finally,  S0  is  evaluated  as  the 
mean  bed  dope  of  the  river  reach. 

The  fundamental  dimensionless  parameters  of  eq 
12  are  the  Courant  number  Cr,  the  Froude  number 
F0,  a  friction  parameter  Fj,  and  the  slope  ratio  S. 
Froude  number  F0  represents  a  ratio  of  speeds,  com¬ 
paring  the  characteristic  flow  velocity  to  the  surface 
wave  celerity  evaluated  at  the  characteristic  flow 


depth.  With  few  exceptions,  the  range  of  Froude 
numbers  for  rivers  is  between  0  and  2,  most  com¬ 
monly  between  0  and  1 .  Slope  ratio  S  compares  the 
river  bed  slope  to  the  energy  gradient  of  the  flow 
evaluated  at  the  characteristic  depth,  velocity,  and 
channel  conveyance.  For  free-flowing  rivers  the 
value  of  the  slope  ratio  is  approximately  1 .  It  is 
higher  for  reservoirs,  where  the  presence  of  back¬ 
water  significantly  reduces  the  energy  gradient. 

The  magnitudes  of  the  terms  in  eq  1 2  directly  re¬ 
flect  the  importance  of  the  physical  processes  they 
represent.  The  behavior  of  solutions  of  eq  12  can  be 
assessed  by  considering  Cr  and  F,  together  with  Fc 
and  Z>],  forming  a  group  of  dimensionless  scaling 
parameters.  Courant  number  Cr  is  composed  of  a 
characteristic  flow  velocity  and  length  and  time 
scales  that  are  characteristic  of  the  wave  motion;  it 
is  the  ratio  of  characteristic  flow  velocity  to  measured 
wave  celerity,  limiting  the  range  of  possible  values  to 
between  0  and  1.  All  terms  in  eq  12  except  b2v*/ 
bt*2  are  functions  of  the  Courant  number.  The 
terms  having  the  Courant  number  as  their  dimen¬ 
sionless  scaling  parameter  follow  from  the  inertia 
terms  of  eq  2.  The  scaling  parameter  representing 
the  largest  inertia  term  cannot  be  larger  than  3  and 
typically  is  about  1. 

The  effects  of  friction  and  bed  slope  on  river  flow 
are  represented  by  the  first-order  terms  with  coeffi¬ 
cients  F,  and  Fc.  The  values  of  these  scaling  param¬ 
eters  relative  to  1  characterize  the  importance  of 
friction  in  the  momentum  balance  and,  therefore, 
define  the  boundaries  and  govern  the  transitions  be¬ 
tween  bulk,  dynamic,  and  gravity  wave  types.  Fc  is 
related  to  the  kinematic  flow  number  obtained  by 
Woolhiser  and  Liggett  (1967)  from  a  direct  normali¬ 
zation  of  eq  2.  Both  measured  wave  celerity  and  di¬ 
mensionless  wavelength  are  components  of  F,  and  Fc, 
though  neither  alone  is  adequate  to  characterize  fric¬ 
tion.  The  magnitudes  of  the  friction  parameters  vary 
linearly  with  dimensionless  wavelength.  The  meas¬ 
ured  wave  celerity  is  a  component  of  the  Courant 
number.  Bulk  flow  cases,  where  the  flow  velocity 
approaches  the  wave  celerity,  correspond  to  a  large 
Courant  number  and  increased  relative  magnitudes 
of  the  friction  parameters.  For  the  opposite  condi¬ 
tion  of  small  flow  velocity  relative  to  measured  wave 
celerity,  the  friction  parameter  magnitudes  decrease 
due  to  their  dependence  on  the  Courant  number,  in¬ 
dicating  a  reduced  role  of  friction  relative  to  other 
processes.  The  other  basic  component  of  F,  and  Fc 
is  the  dimensionless  channel  conveyance,  which  ac¬ 
counts  for  differences  in  the  rate  of  energy  dissipa¬ 
tion  between  different  river  and  flow  situations. 

Fc  is  presented  in  Figure  1  as  a  function  of  the 


Figure  1.  Friction  scaling  parameter  Fc  as  a  / unc¬ 
tion  of  Courant  number  for  a  range  of  dimension¬ 
less  wavelengths  and  a  typical  dimensionless  chan¬ 
nel  conveyance. 


Courant  number  for  a  range  of  dimensionless  wave¬ 
lengths  and  an  assumed  dimensionless  channel  con¬ 
veyance  of  about  9.  Larger  values  of  channel  con¬ 
veyance  displace  the  band  of  curves  to  the  right, 
corresponding  to  smaller  Fc,  with  the  converse  true 
for  smaller  conveyances.  Courant  numbers  greater 
than  0.4  and  dimensionless  wavelengths  greater  than 
4000  are  typical  in  natural  rivers.  For  these  condi¬ 
tions  Fc  is  greater  than  10,  indicating  at  least  an  or- 
der-of-magnitude  dominance  of  friction  over  inertia. 
In  reservoirs,  wavelengths  can  be  short  due  to  hydro- 
power  operations,  and  flow  occurs  at  greater  depths 
and  smaller  velocities  than  in  free-flowing  rivers. 

This  combination  of  conditions  yields  a  small  Cour¬ 
ant  number,  large  channel  conveyance,  and  small  di¬ 
mensionless  wavelength,  and  Fc  for  reservoir  waves 
can  be  significantly  less  than  1 .  When  an  ice  cover 


Figure  2.  Scaling  parameter  D|  as  a  function 
of  Courant  number  for  a  range  of  Froude 
numbers. 


is  present,  the  dimensionless  wavelength  of  a  given 
wave  and  the  relative  magnitudes  of  the  friction 
parameters  increase. 

The  magnitude  of  scaling  parameter  Dv  indicates 
the  importance  of  pressure  gradient  effects  on  the 
wave.  £>|  represents  the  square  of  the  ratio  of  the 
Courant  number  to  the  Froude  number,  and  can  be 
physically  interpreted  as  the  square  of  the  ratio  of 
the  surface  wave  to  the  measured  wave  celerity. 
Figure  2  presents  D,  as  a  function  of  Courant  num¬ 
ber  for  a  range  of  Froude  numbers.  A  value  of  £),  = 
1  Indicates  that  the  measured  wave  celerity  is  equal 
to  the  surface  wave  celerity  at  die  mean  flow  depth. 
When  D,  is  significantly  greater  than  1 ,  measured 
wave  celerity  is  much  smaller  than  the  correspond¬ 
ing  dynamic  wave  celerity,  generally  favoring  bulk 
wave  behavior. 

When  the  friction  parameters  are  much  greater 
than  1 ,  river  flow  waves  exhibit  bulk  wave  behavior 
and,  to  a  good  approximation,  eq  12  becomes 


7*  W+  F  -s)c' 


=  CtDy* 


d2v* 

ax*5-’ 


(13) 


where 


D  = 


£i 


c.1 


2F’(^) 

■'0 


a  nondimensional  form  of  eq  7.  The  diffusion  and 
kinematic  wave  equations  are  contained  within  the 
Saint-Venant  system  equation.  For  high-friction 
flow  conditions,  the  bulk  wave  portion  of  the  system 
equation  is  dominant,  and  the  apparent  wave  celer¬ 
ity  is  that  of  a  bulk  wave.  The  bulk  flow  waves  de¬ 
scribed  by  eq  13  are  independent  of  the  characteris¬ 
tic  equations  3  and  4. 

The  bulk  wave  Courant  number  is  generally  great¬ 
er  than  0.5,  as  the  flow  velocity  is  only  slightly  less 
than  the  wave  celerity.  The  magnitude  of  D  relative 
to  1  indicates  the  importance  of  wave  diffusion  rela¬ 
tive  to  advection.  D  for  bulk  river  waves  is  generally 
less  than  1 ,  indicating  that  the  waves  are  advection- 
dominated.  With  friction  dominating  inertia  and  ad¬ 
vection  dominating  diffusion,  the  equations  describ¬ 
ing  bulk  waves  have  evolved  from  a  second-order  hy¬ 
perbolic  system  to  a  dominantly  hyperbolic  first- 
order  system.  When  S0  the  dimensionless 
diffusion  coefficient  D  varies  inversely  with  the  prod¬ 
uct  of  bed  slope  and  nondimensional  wavelength. 


D  is  presented  in  Figure  3  as  a  function  of  dimension¬ 
less  wavelength  for  a  range  of  Froude  numbers  and  a 
dimensionless  channel  conveyance  of  about  9. 

Dynamic  wave  behavior  is  indicated  when  the  mag¬ 
nitude  of  one  or  both  of  the  friction  parameters  is  on 
the  order  of  1 .  As  friction  and  inertia  both  make  im¬ 
portant  contributions  to  the  momentum  balance,  dy¬ 
namic  wave  behavior  is  an  intermediate  condition  be¬ 
tween  bulk  wave  and  frictionless  gravity  wave  behav¬ 
ior.  Dynamic  wave  behavior  occurs  over  a  wide  range 
of  conditions,  and  wave  subtypes  can  be  identified 
that  correspond  to  varying  magnitudes  of  the  Cour¬ 
ant  number  and  friction  parameters. 

At  the  transition  from  bulk  to  dynamic  waves, 
friction  is  quite  great.  The  magnitudes  of  Fc  and  Cr 
are  both  on  the  order  of  1 ,  and  all  terms  of  eq  1 2  are 
significant,  defining  the  “complete  equation”  case. 

Fc  on  the  order  of  1  and  Cr  «  1  defines  the  “tran¬ 
sition”  case,  where  the  flow  velocity  is  small  relative 
to  observed  wave  speed  and  the  contribution  of 
friction  is  smaller  than  in  the  complete  equation  case. 
The  small  Courant  number  indicates  that  transition 
waves  do  not  retain  significant  bulk  wave  character. 
The  mathematical  description  is  also  simpler,  as  many 
of  the  terms  of  eq  12,  including  the  nonlinear  first- 
order  terms,  are  negligible.  As  friction  is  reduced  fur¬ 
ther,  Fc  and  Cr  are  «  1 ,  and  F,  is  reduced  to  the 
order  of  1 .  Negligible  Fc  indicates  the  lack  of  influ¬ 
ence  of  the  channel  bed  slope  on  the  flow.  These 
conditions  are  representative  of  flow  in  reservoirs. 


Figure  3.  Dimensionless  diffusion  coefficient  of  bulk  waves  as  a  func¬ 
tion  of  dimensionless  wavelength  for  a  range  of  Froude  numbers  and 
a  typical  dimensionless  channel  conveyance. 


The  dimensional  system  equation  for  the  “reservoir” 
case  is 
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where 


7  =  c£r  ’ 

the  familiar  telegraph  equation  or  equation  of  a 
damped  vibration.  Assuming  that  c  and  y  are  con¬ 
stants,  eq  14  has  solutions  of  the  form 


v  =  e-1"  e*!*  *  c  v>  -rVt*1*1)  r] 

where  A:  is  an  arbitrary  nonzero  constant.  The  resis¬ 
tance  term  dv/dr  of  eq  14  follows  from  the  Fl  term 
of  eq  12;  it  provides  exponential  wave  attenuation 
and  retards  wave  celerity,  separating  the  behavior  of 
reservoir  dynamic  waves  from  that  of  gravity  waves. 
The  amplitude  decay  of  reservoir  waves  is  directly 
proportional  to  the  magnitude  of  F,,  indicating  a 
gradual  transition  between  wave  types. 

Gravity  wave  behavior  with  dominant  inertia  oc¬ 
curs  when  friction  is  negligible:  F,  and  Fe  «  1 . 
Gravity  waves,  like  bulk  and  dynamic  waves,  have 
subtypes,  each  of  which  is  undamped  due  to  the 
absence  of  friction.  The  “wave  equation”  case  oc¬ 


curs  when  CT  «  1 .  The  system  equation  for  this 
case  is  the  classical  wave  equation,  obtained  from 
eq  14  by  setting  7  =  0.  Physically,  this  condition  is 
approximated  by  a  smooth,  deep  forebay  of  a  hydro- 
power  station  in  which  the  flow  velocity  is  small  rel¬ 
ative  to  the  surface  wave  speed.  The  “simple  wave” 
case  occurs  when  the  Courant  number  is  on  the  or¬ 
der  of  1 .  The  method  of  characteristics  provides  a 
convenient  formulation  (eq  S,  6)  for  constructing 
solutions  for  this  case.  Positive  disturbances  upstream 
generate  waves  with  converging  characteristics,  lead¬ 
ing  to  wave  steepening  and  shock  formation.  Wave 
propagation  over  small  distances  in  a  smooth  chan¬ 
nel  may  be  adequately  represented  as  a  simple  wave, 
but,  as  shown  in  Figure  1 ,  friction  increases  with 
the  Courant  number  eq  12,  limiting  the  applicability 
of  these  equations. 

The  scaling  parameter  magnitudes  determine 
wave  behavior  and  type.  The  structure  for  under¬ 
standing  river  wave  behavior  developed  from  our 
analysis  is  depicted  in  Figure  4,  where  the  progres¬ 
sion  of  wave  subtypes  is  presented  in  a  natural  order. 
Only  one  dimensionless  parameter  changes  in  mag¬ 
nitude  between  two  adjacent  wave  subtypes,  tend¬ 
ing  from  friction-dominated  kinematic  waves  to 
frictionless  gravity  waves.  The  continuous  depend¬ 
ence  of  the  scaling  parameters  on  the  scaling  variables 
indicates  gradual  transitions  between  wave  types  and 
subtypes. 


Figure  4.  River  wave  structure  obtained  from  the  analysis  of  the  Saint -Venant 
equations. 


Probabilistic  scaling  parameters 

Our  analysis  of  the  equations  governing  river 
waves  provides  a  structure  with  which  to  consider 
wave  behavior.  Choosing  representative  scaling  var¬ 
iables  v0,y0,  Ax,  At,  Cm,  and  S0 ,  we  obtain  quanti¬ 
tative  estimates  of  the  dimensionless  scaling  parame¬ 
ters  F\,FC,  Ct,  and  Dv  For  cases  sufficiently  removed 
from  the  transitions  between  wave  types,  these  in¬ 
tuitive,  deterministic  parameters  may  provide  suffi¬ 
cient  information  to  understand  wave  dynamics. 
However,  assuming  representative  or  mean  scaling 
variables  provides  no  information  about  the  effects 
on  the  wave  of  the  physical  variability  of  the  river. 

It  could  be  argued  that  single-valued  variables  cannot 
meaningfully  represent  natural  rivers,  which  have 
significant  geometric  variability.  An  example  of  this 
problem  is  attempting  to  characterize  reservoir  depth 
by  a  mean  of  6  m,  when  depth  varies  between  2  m 
upstream  and  1 5  m  downstream.  Further  insight 
concerning  the  transitions  between  river  wave  types 
requires  a  more  complete  description  of  the  scaling 
parameters  than  is  provided  by  estimates  of  their 
mean  values.  A  practical  modeling  requirement  that 
has  not  been  addressed  is  how  to  assess  the  confidence 
of  wave  type  interpretation  when  scaling  parameters 
fall  near  a  transition. 

We  will  now  consider  the  more  general  case  where 
scaling  variables  are  continuous,  possessing  a  mean 
and  a  distribution.  By  treating  these  physical  quan¬ 
tities  as  random  variables,  probabilistic  concepts  can 
be  applied.  As  functions  of  random  variables,  the 
scaling  parameters  are  themselves  random  variables. 
Rosenblueth  (1975)  developed  a  method  for  estimat¬ 
ing  the  first  and  second  moments  and  minimum  val¬ 
ues  of  a  function  of  random  variables,  given  the  first 
two  or  three  moments  of  each  variable  and  correla¬ 
tion  coefficients  for  each  pair  of  random  variables. 

The  estimate  obtained  of  the  mean  value  of  the  func¬ 
tion  is  equivalent  to  a  second-order  Taylor  series  ap¬ 
proximation,  and  the  estimate  of  the  coefficient  of 
variation  is  a  first-order  approximation.  The  method 
is  algebraic,  replacing  the  distribution  of  each  random 
variable  by  point  estimates  and  not  requiring  the  com¬ 
putation  of  derivatives. 

Before  proceeding,  the  Rosenblueth  method  for 
a  function  of  two  correlated  random  variables  will 
be  described  in  detail.  The  coefficient  of  variation 
of  a  random  variable  x  is  defined  as 

(16) 

where  Sx  is  the  standard  deviation  of  x  and  x  is  the 
mean.  The  correlation  coefficient  of  random  vari¬ 
ables  xt  and  x2  is  defined  as 


.=cov(x1£a) 
p  S,  S,  • 


(17) 


where  cov(xt  ,x2 )  is  the  covariance.  A  correlation 
coefficient  of  magnitude  1  indicates  a  perfect  linear 
correlation  between  the  variables,  while  a  coefficient 
of  0  indicates  uncorrelated  variables.  The  scaling 
variable  ranges  corresponding  to  the  propagation  of 
a  river  wave  through  a  reach  are  relatively  small. 
Numerical  experiments  have  shown  that  if  the  rela¬ 
tionship  between  random  variables  over  a  limited 
range  can  be  written  as 

x,  =ax5,  (18) 

then  the  magnitude  of  p  will  be  approximately  1 . 

In  general,  sufficient  data  will  not  be  available  to 
describe  the  precise  nature  of  the  distribution  of 
each  scaling  variable.  Therefore,  for  simplicity  we 
will  assume  that  these  distributions  are  symmetric 
about  the  mean  value.  Hie  point  estimates  repre¬ 
senting  the  distribution  of  a  scaling  variable  are  then 


(19) 


These  point  estimates  sum  to  1  and,  in  the  case  of 
uncorrelated  random  variables,  are  each  Va.  The 
point  estimates  provide  weighting  factors  for  evalu¬ 
ating  the  function  at  points  a  standard  deviation 
from  the  mean  of  each  random  variable: 


y++=y(x i  +sXj  ,xa  +sXj), 

y+-=y(x i  +«**,  ~sXt), 

y.+  =;K*i  ~sXl  +sXt), 

y~.=y(x i  -sXj  ,x2  -sXt). 


(20) 


Then,  the  expected  values  or  moments  of  the  func¬ 
tion  can  be  obtained: 


£lvn]  =  F++y;++?+^_ 
+  P_+yn_+  +F__y2_  . 


(21) 


The  expected  value  or  mean  of  y  is  found  by  setting 
n  =  1.  The  variance  ofy  is  readily  computed  as 


Ktyl-fly2]-^])1  .  (22) 
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Table  1.  Application  of  Rosenblueth  (1975)  method  to 
Liard  River  data. 


Input:  v0=  0.72  m /»,  KV(  =  0.25,  c  =  1.4  m/i,  Vc  =  0.20 
Output:  y  =  Ct  =  O.S3S7(++),  0.8036(+-),  0.3214(-+),  0.482 1(— ) 


0  =  0,6 

0  =  0.8 

0  =  1.0 

*lcr] 

0.S20 

0.514 

0.S09 

Hcr] 

0.0127 

0.0067 

0.0007 

0.216 

0.159 

0.053 

r 


Rosenblueth  (1975)  discussed  the  generalization  of 
this  technique  to  functions  of  any  number  of  ran¬ 
dom  variables. 

A  sample  calculation  using  Rosenblueth ’s  proce¬ 
dure  to  obtain  parameters  describing  the  Courant 
number  distribution  for  data  from  the  Liard  River 
(Parkinson  1981 , 1982)  is  given  in  Table  1 .  The 
coefficients  of  variation  chosen  for  wave  celerity 
bound  most  of  the  observed  values  within  one  stan¬ 
dard  deviation  of  the  mean.  Because  the  mean  flow 
velocity  is  a  computed  value,  a  relatively  large  coef¬ 
ficient  of  variation  for  velocity  was  assumed.  Row 
velocity  and  wave  celerity  for  bulk  flow  waves  can 
be  related  by  an  equation  of  the  form  of  eq  1 8  and 
p  *  I.  However,  to  examine  the  sensitivity  of  the 
computation  to  the  correlation  coefficient,  results 
are  given  in  Table  1  for  three  different  p  values. 

With  increasing  p  the  mean  Courant  number  decreas¬ 
es  only  slightly,  but  the  variance  and  coefficient  of 
variation  decrease  significantly.  Estimates  of  the 
maximum  and  minimum  values  of  the  Courant  num¬ 
ber  distribution  should  be  greater  than  and  less 
than  Cr  + ,  respectively. 

Given  a  function  of  random  variables,  Rosen- 
blueth’s  method  provides  an  improved  estimate  of 
the  mean  and  estimates  of  scatter  about  the  mean 
and  the  limits  of  the  distribution.  The  mean,  stan¬ 
dard  deviation,  and  limits  of  a  random  variable  unique¬ 
ly  specify  a  0  distribution.  Harr  (1977)  gives  the 
general  expression  for  this  distribution  as 

^)  =  (h-«)B(<H-l,0+l)(fcfl)  fea)0  (23) 

where 


and 


where  a  and  b  are  the  minimum  and  maximum  val¬ 
ues  of  y,  respectively,  and  B  is  a  beta  function. 

When  a  and  0have  the  same  sign,  the  0  distribution 
is  unimodal  and  bell-shaped.  The  coefficients  of 
skewness  and  kurtosis  of  the  0  distribution  can  be 
readily  obtained.  With  this  distribution  the  proba¬ 
bility  of y  in  a  given  range  can  be  determined. 

To  incorporate  the  effects  of  the  physical  varia¬ 
bility  of  the  river  on  the  wave,  we  will  use  the  Ros¬ 
enblueth  method  and  obtain  a  corresponding  0  dis¬ 
tribution  of  the  friction  scaling  parameters.  This 
procedure  allows  further  investigation  of  the  bulk/ 
dynamic  and  dynamic/gravity  wave  transitions  and 
provides  a  reliability  estimate  for  wave  type  deter¬ 
mination  when  the  friction  parameter  mean  is  near 
a  transition. 

Case  studies 

The  river  wave  structure  that  has  been  developed 
depends  on  parameter  magnitudes  that  are  large  or 
small  relative  to  1  or  are  on  the  order  of  1 .  We  will 
now  use  experience  from  case  studies  to  define  more 
precisely  the  parameter  ranges  corresponding  to  each 
wave  type  and  transition.  We  will  also  classify  a  num¬ 
ber  of  diverse  river  wave  cases  and  consider  changing 
wave  behavior  over  distance  to  test  the  usefulness 
and  practicality  of  our  analysis. 

Individual  case  studies  and  physical  variables  rep¬ 
resentative  of  each  river  and  wave  are  presented  in 
Table  2,  as  well  as  the  first-order  Taylor  series  ap- 


«  =  (l-yMltv)  , 


0  =  ^--(a+2)  , 


proximations  of  the  mean  values  of  the  scaling  param- 
eters-that  is,  the  parameter  evaluated  at  the  mean 
values  of  all  scaling  variables.  Estimated  Manning’s 
roughness  and  river  bed  slope  values  are  designated 
in  the  table  with  an  e.  For  free-flowing  rivers  where 
measured  velocities  were  unavailable,  die  Chezy 
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equation  with  the  characteristic  depth,  river  bed 
slope,  and  conveyance  was  used  to  obtain  a  comput¬ 
ed  mean  characteristic  velocity,  which  is  designated 
with  a  c.  The  mean  velocity  for  reservoirs  was  com¬ 
puted  using  the  continuity  equation  from  the  mean 
discharge,  channel  width,  and  characteristic  depth. 

Bulk  flow  waves  occur  when  the  friction  parame¬ 
ters  are  sufficiently  large  relative  to  1 .  Because  the 
Courant  number  is  by  definition  less  than  1 ,  Fc  is 
always  smaller  than  Fx  and  is  the  natural  parameter 
to  define  the  boundary  between  bulk  and  dynamic 
waves.  Parkinson  (1981 , 1982)  studies  the  river 
waves  that  initiated  ice-cover  breakup  on  the  Liard 
and  Mackenzie  Rivers.  These  waves  had  long  peri¬ 
ods,  characteristic  of  spring  flood  waves  in  large 
river  systems.  The  bulk  wave  behavior  that  occurred 
is  reflected  in  the  large  friction  parameter  magnitudes 
for  these  cases.  Large  values  of  Fc  also  correspond 
to  data  from  open-water  and  flow-over-the-ice  con¬ 
ditions  in  a  free-flowing  reach  of  the  Ottauquechee 
River  (Ferrick  and  Lemieux  1983),  again  indicating 
bulk  wave  behavior.  Ferrick  (1980)  numerically 
modeled  rapidly  varied  flow  waves  from  laboratory 
studies  of  dam  breaches  at  the  Waterways  Experi¬ 
ment  Station.  A  wave  release  to  a  previously  dry 
channel,  corresponding  to  a  mean  Fc  value  of  1 1 , 
was  precisely  simulated  using  a  bulk  wave  numerical 
model.  In  addition,  Ferrick  et  al.  (1984)  found  that 
rapidly  varied  flow  waves  in  the  Hiwassee  and  Clinch 
Rivers,  with  a  minimum  mean  Fc  value  of  12,  behaved 
as  bulk  waves. 

According  to  these  observations,  Fc  greater  than 
about  10  indicates  bulk  wave  behavior.  The  system 
equation  terms  that  are  more  than  an  order  of  mag¬ 
nitude  smaller  than  the  other  terms  are  apparently 
negligible.  With  this  guideline  we  can  approximate 
the  transitions  between  wave  types  and  subtypes  as 
corresponding  to  scaling  parameter  magnitudes 
greater  than  10  (»  1)  and  less  than  0.1  («  1).  Be¬ 
cause  the  scaling  parameters  representing  a  given 
river  reach  and  wave  are  continuous,  however,  any 
single-valued  definition  of  a  transition  is  inadequate. 
The  governing  equations  used  for  analysis  change 
between  bulk,  dynamic,  and  gravity  wave  conditions, 
and  these  transitions  will  be  considered  in  detail. 

All  of  the  bulk  wave  cases  have  mean  Courant 
numbers  greater  than  0.S  (Table  2),  a  necessary  but 
apparently  not  sufficient  condition  for  the  occur¬ 
rence  of  bulk  waves.  Table  2  also  confirms  that  large 
dimensionless  wavelengths  favor  but  do  not  guaran¬ 
tee  bulk  wave  behavior.  The  wave  diffusion  present 
in  all  cases  was  sufficient  to  prevent  the  formation 
of  a  shock  wave  or  depth  discontinuity  at  the  wave 
front.  However,  as  the  magnitude  of  D  approaches 
0,  the  probability  of  wave  steepening  and  shock 


formation  increases.  Only  in  the  Clinch  River  was 
wave  diffusion  sufficient  to  produce  significant 
wave  peak  attenuation.  These  observations  support 
the  order-of-magnitude  approximation  for  locating 
the  transition  between  the  bulk  wave  subtypes.  Mean 
dimensionless  diffusion  coefficients  less  than  about 
0.1  indicate  kinematic  waves,  and  diffusion  waves 
occur  when  D  is  larger. 

The  case  studies  from  Table  2  that  bracket  the 
bulk/dynamic  wave  transition  Fc  «  10  are  waves 
having  a  1-hr  period  on  the  Clinch  River  and  a  short- 
period  wave  on  the  Ottauquechee  River  that  occurred 
immediately  following  ice  cover  breakup  (moving 
ice).  In  this  discussion  we  refer  to  these  cases  as 
simply  Clinch  River  and  Ottauquechee  River.  The 
extensive  documentation  available  from  the  Clinch 
River  study  allows  well-supported  estimates  of  the 
coefficient  of  variation  for  the  scaling  variables. 

These  values  were  chosen  so  that  the  majority  of  the 
data  fall  within  one  standard  deviation  of  the  mean, 
and  essentially  ail  data  are  within  three  standard 
deviations  of  the  mean.  With  less  data  available 
from  the  Ottauquechee  River,  we  assume  for  pur¬ 
poses  of  comparison  that  the  coefficients  of  varia¬ 
tion  equal  those  from  the  Clinch  River.  The  coeffi¬ 
cients  of  variation  of  velocity,  wave  celerity,  wave¬ 
length,  and  hydraulic  depth  were  taken  as  0. 1 5 ,  and 
as  0.20  for  channel  conveyance.  For  both  cases, 
correlation  coefficients  for  velocity  /depth,  convey¬ 
ance/depth,  velocity/celerity,  and  wavelength/celer¬ 
ity  were  specified.  Observed  bulk  wave  celerity  on 
the  Clinch  River  indicates  a  velocity /wavelength 
correlation.  The  relatively  high  wave  celerity  on  the 
Ottauquechee  River  indicates  dynamic  wave  effects, 
and  wave  celerity/depth  and  wavelength/depth  cor¬ 
relations.  The  correlated  variables  for  river  flow  are 
related  analytically  in  the  form  of  eq  18,  and  there¬ 
fore  correlation  coefficients  of  1  were  assumed. 

With  these  input  data,  Rosenblueth’s  method 
yields  the  mean,  coefficient  of  variation,  and  mini¬ 
mum  and  maximum  Fe  values  presented  in  Table  3. 
Table  3  also  includes  parameters  of  the  ^-distribution 
fit  to  Fc  and  probabilities  of  Fc  in  selected  ranges  for 
each  case.  We  note  that  the  second-order  estimates 
of  the  mean  for  each  case  are  significantly  larger 
than  the  first-order  estimates  in  Table  2.  The  rela¬ 
tively  large  coefficients  of  variation  of  Fc  indicate 
that  these  distributions  have  significant  spread.  The 
0  distributions  representing  the  Fe  distributions  for 
the  Clinch  and  Ottauquechee  Rivers  (Fig.  5)  extend 
a  significant  distance  on  either  side  of  the  approxi¬ 
mate  transition,/^  =  10.  Therefore,  a  finite  proba¬ 
bility  exists  that  bulk  or  dynamic  wave  behavior  will 
occur  in  each  case.  This  fact  implies  that  point  esti¬ 
mates  of  Fc  falling  near  a  transition  are  of  limited 
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Figure  5.  Probability  distributions  of  friction  scaling  parameter 


Fc  representing  case  studies  of  the  Clinch  River  and  the  Ottau- 


quechee  River. 


value  as  indicators  of  wave  type,  and  emphasizes 
that  the  transitions  between  wave  types  in  rivers  are 
continuous.  If  we  consider  the  transition  range 
8  <FC  <  12, the  ratio  of  probabilities,  Rt ,  of  bulk- 
to-dynamic  wave  behavior  for  Fc  outside  this  range  is 


= 


P[FC  >  12] 


P[FC  <  8] 


(24) 


Values  of  /?i  near  1  do  not  identify  wave  type,  but 
progressively  larger  deviations  from  1  indicate  wave 
type  with  increasing  reliability.  In  our  case,/?, 
(Clinch)  *  3.6  and  Rt  (Ottauquechee)  =  0.35,  indi¬ 
cating  probable  bulk  wave  behavior  and  probable  dy¬ 
namic  wave  effects,  respectively.  The  diagnosis  of 
Clinch  River  bulk  wave  behavior  is  strongly  support¬ 
ed  by  modeling  results  (Ferrick  et  al.  1 984). 

Dynamic  waves  span  a  range  of  the  friction  param¬ 
eters  of  about  two  orders  of  magnitude.  The  scaling 
parameter  bounds  that  identify  complete-equation 


dynamic  waves  are  roughly  defined  as  Fe  <  10  and 
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Q  >  0.1.  In  these  relatively  high  friction  cases, 
where  dynamic  waves  typically  result  from  some 
rapidly  varying  initial  condition,  frictional  effects 
can  eventually  restore  bulk  wave  behavior.  Down- 
stream-propagating  waves  following  a  dam  breach 
or  the  sudden  release  of  an  ice  jam  are  examples  of 
such  cases.  Stoker  (1957)  simulated  the  propagation 
of  the  front  of  a  hypothetical,  rapidly  varying,  large- 
amplitude  flood  wave  on  the  Ohio  River  using  the 
dynamic  wave  equations.  At  the  upstream  boundary 
he  imposed  a  1 ,5-m/hr  rise  in  the  water  level  for  a 
4-hr  period  between  initial  and  final  steady-flow 
conditions.  The  magnitudes  of  the  friction  parame¬ 
ters  and  the  Courant  number  for  2  hr  after  the  initi¬ 
ation  of  the  wave  (Table  2)  indicate  significant  iner¬ 
tia.  However,  average  scaling  variables  representative 
of  the  first  10  hr  of  the  simulation  yield  mean  fric¬ 
tion  parameters  in  the  bulk  wave  range.  Mean  scal¬ 
ing  variables  representing  the  latter  part  of  this  time 
period  produce  even  larger  friction  parameter  magni¬ 
tudes  and  a  mean  Courant  number  greater  than  0.5. 


Stoker  stated  that  prior  to  his  work  more  gradual 
flood  waves  recorded  on  the  Ohio  River  had  been 
successfully  described  with  bulk  wave  models,  as 
would  be  expected  based  on  our  analysis. 

The  other  complete-equation  dynamic  wave 
cases  given  in  Table  2  are  predominantly  rapidly 
varying  flow  waves  with  small  dimensionless  wave¬ 
lengths.  This  group  of  cases  includes  a  laboratory 
dam  breach  release  to  a  channel  with  a  significant 
base  flow  (WES  flume-wet),  laboratory  studies  (Fer- 
rick  1984)  of  the  effects  of  bed  slope  and  ice  cover 
upon  a  sequence  of  flow  waves  (CRREL  flume), 
and  three  waves  on  the  Ottauquechee  River  immedi¬ 
ately  prior  to  and  during  ice  cover  breakup.  The 
laboratory  flume  cases  each  had  bulk  measured  wave 
speeds,  but  because  of  very  smooth  channels  the 
friction  parameters  fall  in  the  dynamic  wave  range. 
At  break-up  the  Ottauquechee  study  reach  was 
ponded  to  a  large  depth  relative  to  wave  amplitude. 
The  channel  conveyance  at  the  breakup  was  much 
smaller  than  in  the  laboratory  studies,  but  the  meas¬ 
ured  river  wave  celerities  approached  those  of  dy¬ 
namic  waves  and  are  reflected  in  small  Courant  num¬ 
bers.  The  upstream  half  of  ice-covered  Pine  Falls 
reservoir  (Kartha  1976)  is  a  complete-equation  dy¬ 
namic  wave  case  that  did  not  involve  a  short-period 
flow  wave.  However,  a  large  channel  conveyance 
and  a  low  flow  velocity  relative  to  measured  wave 
celerity  compensate  for  a  relatively  large  dimension¬ 
less  wavelength. 

The  Courant  number  for  the  transition  dynamic 
wave  case  is  smaller  than  that  for  the  complete  equa¬ 
tion  case,  and  the  friction  parameter  Fc  is  also  gener¬ 
ally  smaller.  The  dimensionless  scaling  parameter 
bounds  that  identify  transition  dynamic  waves  are 
roughly  defined  as  Fc>  0.1  and  Cr  <  0.1 .  Waves  of 
2-hr  and  6-hr  durations  in  the  ice-covered  St.  Lawr¬ 
ence  River  between  Lake  Ontario  and  Cardinal, 
Ontario  (Yapa  1983),  are  examples  of  transition-case 
dynamic  waves.  The  observed  wave  celerity  was  that 
of  a  dynamic  wave  and  was  much  greater  than  the 
flow  velocity. 

The  reservoir  dynamic  wave  case  is  characterized 
by  a  small  flow  velocity  relative  to  measured  wave 
speed  and  relatively  small  frictional  effects.  This 
case  is  roughly  defined  by  dimensionless  parameters 
Fc  and  Ct  <0.1.  In  Table  2  waves  from  hydroelec¬ 
tric  power  generation  in  Old  Hickory  Reservoir  on 
the  Cumberland  River  (Ferrick  1979),  in  the  down¬ 
stream  half  of  Pine  Falls  Reservoir  on  the  Winnipeg 
River  (Kartha  1976),  and  in  Wheeler  Reservoir  on 
the  Tennessee  River  (Ferrick  and  Waldrop  1977)  are 
reservoir  waves.  Wave  propagation  in  the  upstream 
portion  of  the  Old  Hickory  and  Wheeler  Reservoirs 
is  subject  to  larger  frictional  effects  than  farther 


downstream  and,  if  analyzed  separately,  would  be 
classified  as  a  higher-friction  dynamic  wave  type. 

For  each  of  these  cases  friction  remains  important, 
as  velocity  decay  to  1 0%  of  the  original  amplitude 
occurs  within  a  4-  to  6-hr  period  (eq  1 S). 

The  case  studies  compiled  in  Table  2  do  not  in¬ 
clude  a  case  that  falls  in  the  gravity  wave  range.  The 
reservoir  dynamic  waves  are  the  smallest  friction 
cases  available,  with  estimated  mean  Ft  values  of 
about  1  from  the  Old  Hickory,  Pine  Falls,  and 
Wheeler  Reservoir  data.  We  will  use  the  2-hr  wave 
period  data  from  Wheeler  Reservoir  to  consider  the 
dynamic/gravity  wave  transition.  A  treatment  of 
this  boundary  consistent  with  that  used  for  the 
bulk/dynamic  wave  transition  provides  values  of  Fj 
between  0.083  and  0.1 25  as  the  transition  range. 
Table  4  presents  wave  amplitude  decay  times  of 
2-hr  period  waves  on  Wheeler  Reservoir  for  Ft  as 
computed  and  bounding  the  transition.  The  actual 
conditions  are  far  removed  from  the  transition. 

The  scaling  variables  presented  in  Table  2  for 
Wheeler  Reservoir  are  mean  values  for  the  complete 
reservoir.  These  values,  repeated  in  Table  S  along 
with  estimated  coefficients  of  variation,  provide  in¬ 
put  for  the  Rosenblueth  method.  Wheeler  is  a  typi¬ 
cal  reservoir  in  that  dramatic  changes  in  depth  and 
flow  velocity  occur  over  its  length,  reflected  in  large 
coefficients  of  variation.  If  we  partition  the  reser¬ 
voir  and  consider  the  40-km  reach  farthest  down¬ 
stream,  this  physical  variability  is  greatly  reduced. 
AD  coefficients  of  variation  are  reduced  by  roughly 
a  factor  of  2,  which,  for  purposes  of  comparison. 

Table  4.  Wave  amplitude  decay  time  (hr). 
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Table  5.  Input  data  for  Rosenblueth ’s  method. 
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Table  6.  Dynamic/gravity  wave  transition. 
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Figure  6.  Probability  distributions  of  friction  scaling  parame¬ 
ter  F /  representing  case  studies  of  Wheeler  Reservoir  and  the 
downstream  portion  of  Wheeler  Reservoir. 


we  assume  is  exactly  the  case.  In  addition,  the  mean 
flow  depth  in  the  lower  portion  of  the  reservoir  is 
greater  and  the  velocity  is  smaller  than  values  repre¬ 
senting  the  complete  reservoir,  moving  Fj  toward 
the  dynamic/gravity  wave  transition.  These  data  are 
also  given  in  Table  5.  The  correlation  coefficients 
used  in  the  method  were  identical  to  those  for  dy¬ 
namic  waves  in  the  Ottauquechee  River. 

Presented  in  Table  6  are  parameters  of  the  Fj 
distribution  for  both  cases  from  the  Rosenblueth 
method,  parameters  describing  a  ^-distribution  fit 
to  Fj,  and  probabilities  bounding  and  within  the  dy¬ 
namic/gravity  wave  transition.  The  second-order 
estimate  of  the  mean  for  the  complete  reservoir  is  a 
third  smaller  than  the  first-order  estimate  given  in 
Table  2.  The  mean  F|  of  the  downstream  portion 
is  about  one-half  the  value  for  the  complete  reser¬ 
voir,  with  a  significantly  smaller  coefficient  of  varia¬ 
tion.  The  ratio 

F[F,>  0.125]  ,  % 

R*  ~F[F,  <0.083]  ’  (25) 


where  a  large  value  of  R2  relative  to  1  indicates  dy¬ 
namic  wave  behavior  and  an  R2  value  much  less 
than  1  indicates  gravity  wave  behavior,  reveals  that 
neither  case  approaches  the  wave  transition.  The 
pronounced  difference  in  spread  between  the  F\ 
distributions  and  the  small  probability  that  FI  is 
less  than  0.1  can  be  readily  seen  with  a  plot  of  these 
distributions  (Fig.  6).  These  results  demonstrate 
that  river  waves  in  reservoirs  may  resemble  gravity 
waves  over  small  distances,  but  waves  propagating 
more  than  a  few  wavelengths  are  significantly  af¬ 
fected  by  friction. 


CONCLUSIONS 

The  Saint-Venant  system  equation,  formed  by 
combining  the  continuity  and  momentum  equations, 
contains  a  pair  of  wave  celerities.  The  dynamic 
wave  celerity  is  associated  with  the  second-order 
terms,  which  are  due  to  inertia  and  the  pressure 
gradient,  and  the  kinematic  wave  celerity  is  associ- 


ated  with  the  first-order  terms,  which  are  due  to  bed 
slope  and  friction.  River  wave  behavior  is  controlled 
by  the  balance  between  friction  and  inertia.  The 
nondimensional  scaling  parameters  of  the  Saint- 
Venant  system  equation  can  be  used  to  quantify 
this  balance.  The  dimensionless  parameters  define 
a  spectrum  of  river  waves,  with  continuous  transi¬ 
tions  between  wave  types  and  subtypes.  Dynamic, 
gravity,  diffusion,  and  kinematic  waves  correspond 
to  specific  scaling  parameter  ranges  of  this  spectrum. 
Bulk  waves  occur  when  the  friction  terms  of  the 
system  equation  dominate  the  inertia  terms.  Dynam¬ 
ic  waves  occur  when  the  second-order  inertia  terms 
of  the  system  equation  are  significant.  Because  fric¬ 
tion  in  rivers  is  generally  important,  there  are  few 
cases  where  the  inertia  terms  dominate  the  friction 
terms  and  river  waves  behave  as  gravity  waves. 

The  parameter  range  corresponding  to  each  wave 
type  and  transition  was  defined  using  data  from 
case  studies,  allowing  direct  application  of  the  scal¬ 
ing  parameters  for  identification  of  wave  type.  The 
capability  to  identify  wave  type  is  necessary  to  the 
construction  of  appropriate  mathematical  models  of 
river  flow.  Changes  in  wave  behavior  with  propaga¬ 
tion  distance  for  rapidly  varying  initial  flow  condi¬ 
tions  and  changes  resulting  from  the  presence  of  an 
ice  cover  can  be  addressed  quantitatively  using  the 
scaling  parameters.  Interpreting  the  scaling  parame¬ 
ters  as  random  variables  supplies  the  generality  need¬ 
ed  to  consider  waves  in  rivers  with  significant  physi¬ 
cal  property  ranges  and  provides  a  measure  of  the 
reliability  of  wave  type  assessments  near  a  transition. 
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